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SUMMARY 

Expressions are presented for the load distribution on a represent- 
ative group of plan forms in sideslip at supersonic speeds. These 
expressions were obtained by the application of lifting-surface theories 
based on the linearized equation for compressible flow. Sketches of the 
load distributions are included. 


INTRODUCTION 

In three recent reports (references 1, 2, and 3 ) the variations of 
rolling moment, of yawing moment, and of lift and pitching moment with 
sideslip have been investigated for a group of wing plan forms for super- 
sonic speeds. The pressure distributions required to compute these 
forces and moments were calculated using linearized compressible— flow 
theory for thin airfoils. Since the reports referred to were concerned 
with the detailed expressions of moments and forces for the various plan 
forms, it was decided that the reference value of the pressure distribu- 
tions, their possible utility in stress analysis and design, and the 
desirability of including some pictorial representations justified the 
treatment of these distributions as the subject of a separate report. 

By virtue of the many approximations involved in its derivation, 
the linearized theory applied constitutes one of the most simplified 
analytical approaches to compressible— flow problems. Furthermore, in 
addition to the* factors approximated in the linearization of the poten- 
tial theory, the analysis employed does not account for the lack of 
complete rigidity of a wing nor the effects of viscosity in the flow. 
These are two important factors that may have considerable effect on 
the actual distribution of the pressure on a wing. Thus, it is not 
expected that these pressure distributions will conform precisely to 
those obtained in the actual physical flow. It is expected, however, 
that these theoretical solutions will be good first approxinlations for 
the plan forms and conditions considered herein and they should provide 
satisfactory indications of the pressure— difference contours in general 
if not in detail. 


2 


NACA TN 2007 


The investigation covers the following configurations: (See figs. 

1 and 2 0 ) (l) Triangular plan forms with subsonic leading edges or with 

supersonic leading edges; (2) trapezoidal plan forms with all possible 
combinations of raked— in, raked— out, subsonic, or supersonic tips; (3) 
rectangular plan forms; and (^-) two swept— back plan forms with super- 
sonic trailing edges developed from the triangular wings. Illustrations 
are included in order to provide a convenient visual correlation between 
the expressions for the pressure distributions and for the moments and 
forces that were calculated from them (references 1, 2, and 3) • The 
arrangement of the appendixes was based on the desire to present a syste- 
matic and convenient compilation of expressions and illustrations for the 
load distributions for the various plan forms considered. 
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SYMBOLS, 


aspect ratio 



COEFFICIENTS, AND AXES 


span of wing measured normal to plane of symmetry 
Mach number parameter 


./M^-l) 


ratio of tangent of right tip angle to tangent of Mach 

cone angle ( — 2 — 

\ tan p 

chord of wing in plane of symme try- 

incomplete elliptic integral of the second kind with 
modulus k ^ J' 1— k 2 sin 2 0 d0 ^ 
complete elliptic integral of the second kind with modulus k 
E (f,k) 

incomplete elliptic integral of the first kind with modulus k 
'<P d0 


) 


___ J 1— k 2 sin 2 0 

complete elliptic integral of the first kind with modulus k 
F (J,k) 


over-all longitudinal length of swept-back wing 

slope of right wing tip measured in plane of wing (positive 
for raked— out tip, negative for raked— in tip) 
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free— stream Mach number 

slopes of plan— form edges relative to wind axes 
pressure differential across wing surface , positive upward 
loading coefficient 


free— stream dynamic pressure 
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area of wing 

perturbation velocity parallel to positive x axis 
free— stream velocity 

perturbation velocity parallel to z axis (positive upward) 
rectangular coordinates of wind axes (fig. 3) 
angle of attack, radians 

sideslip angle (positive when sideslipping to right), degrees 

rectangular coordinates of stability axes (fig. 3) 
rectangular coordinates of body axes (fig. 3) 
air density in the free stream 
perturbation velocity potential 


Mach angle 


( tan 1 




Subscripts 

A,B, . . .V expressions given in Appendix B 

The body axes are generally a right -handed system of three orthog- 
onal axes as shown in figure 3 with the longitudinal axis !' lying in 
the plane of the wing. The stability axes are, in effect, the body axes 
rotated about the lateral axis ri * (through -a.) until the longitudinal 
axis is in the horizontal plane containing the free— stream vector; a 
subsequent rotation about the vertical axis £ (through 0) would bring 
the longitudinal axis in line with the free— stream vector and the axes 
would now be coincident with the wind axes. 
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It should he mentioned that the orientation of these axes as shown 
in figure 3 is convenient for the calculations required to determine 
the pressure distributions and the resulting forces and moments. In the 
application of these results to the calculation of the motion and dynamic 
stability of an airplane, however, the axes are usually rotated so that 
the positive direction of the longitudinal axis is into the free stream 
and the positive direction of the vertical axis is downward, that is, 
toward the undercarriage of the airplane . 


The development of the expressions for the pressure (i»e. ; load) 
distributions on wing plan forms in sideslip was merely an application 
of supersonic wing theory. In this report, only the part of the theory 
relating to the flat— plate or so-called Additional" loading, which is 
the loading resulting from a change in angle of attack, will be consid- 
ered. The loadings due to camber and twist are usually assessed inde- 
pendently, and the sum of these two loadings is often referred to as the 
"basic" load distribution. 

By linearization of the partial differential equation for compres- 
sible flow it is possible to develop a simplified lifting— surface theory 
for thin airfoils. The linearization is made possible by the assumption 
that, for thin airfoils, the perturbation velocities induced by the air- 
foil are small relative to the free— stream velocity. If the free— stream 
velocity vector is parallel to and in the direction of the positive x 
axis and if cp denotes the perturbation velocity potential for i sen- 
tropic flow, the linearized partial differential equation for steady- 
state conditions at supersonic velocities is 


where Mi is the Mach number of the free stream. There have been a 
number of methods developed that provide means of fitting solutions of 
this equation to the boundary conditions of thin— airfoil theory (e.g«, 
references 4 through 9). The results to be given herein were determined 
through the general use of source— sink and doublet distributions (refer- 
ences 4, 5 , 6, and 9)* In particular, the method of reference 6 was 
applied to cases where a subsonic tip occurs in conjunction with a 
supersonic leading edge or tip; whereas the load distributions for all 
other edge and Mach cone arrangements were calculated by application of 
the methods summarized in reference 9- 

The first step in the analysis is the establishment of the boundary 
conditions. For thin airfoils the boundary conditions are usually 
restricted to the z=0 plane. Thus, if the local angles of attack- at 
various spanwise stations of the wing are specified and it is assumed 


METHOD OF ANALYSIS 
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that the wing is coincident with the z=0 plane , the boundary conditions 
are set. Next, the expressions for the loading — and the angle of 

attack a are formulated in terms of parameters that can be related to 
the potential solutions of the differential equation and to the boundary 
conditions. For linearized theory these relationships are 


AP 
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Thus the problem is reduced to determining cp in such a manner that 
— is equal to the specified local angle of attack at every spanwise 
station of the wing. 


The general problem of specifying the angle of attack and of solv- 
ing for the resulting velocity potential is one that usually requires 
the solution of an integral equation. (See reference 9.) For cases 
where the edges of the plan form are supersonic, however, the lack of 
interaction between the upper and lower- surfaces of the wing permits the 
problem to be solved by a distribution of sources in accordance with the 
local slopes of the plan form and a straightforward integration of their 
potentials. The triangular and trapezoidal plan forms with supersonic 
edges were treated in this manner. Likewise, wherever a subsonic edge 
is in conjunction with a supersonic leading edge or tip a straightforward 
integration can be employed. For this case reference 6 provides a method, 
based on the consideration of the upwash between the subsonic edge and the 
Mach cone, whereby the usual operations involved in the solution of the 
integral equation are eliminated. In general, however, it is necessary 
to go through rather involved procedures to calculate the load distribu- 
tion when the camber, twist, and angle of attack of the plan form are 
specified. These procedures are discussed in reference 9 wherein, for 
conical— flow conditions, a loading element is used to set up the inte- 
gral equation and then the usual integral— equation techniques* are 
employed to solve it. 


It should be mentioned also that, when the trailing edge is sub- 
sonic, an additional stipulation based on some physical concept for the 
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flow, such as the Kutta condition which is applied herein, is needed in 
order to eliminate all but one of an infinite number of potential solu- 
tions that will satisfy the boundary conditions; 


PRESENTATION OP RESULTS 

All the plan forms and conditions investigated in references 1, 2, 
and 3 were made up using the five following combinations of straight 
edges: 

1. A supersonic leading edge in conjunction with a subsonic 

leading edge 

2. A supersonic leading edge in conjunction with a subsonic 

trailing edge 

3. Two subsonic leading edges 

4. A subsonic leading edge in conjunction with a subsonic 

trailing edge 

5. Two supersonic leading edges 

The expressions for the load distributions on these five combina- 
tions are given in Appendix A in terms of the wind-axes notation. 

In order to provide an easy correlation between the load distribu- 
tions and the aerodynamic characteristics of the plan forms presented 
in references 1, 2, and 3> Appendix B contains the expressions for the 
load distributions on various sectors of the plan forms in terms of the 
plan— form parameters and the body-axes notation. Since the plan forms 
are restricted to the z=0 plane for the purpose of analysis by the 
thin-airfoil theory, the body-axes notation, in a sense, refers to coor- 
dinates on the projection of the plan form onto the z=0 plane, which 
corresponds to giving the coordinates in terms of the stability-axes 
notation. This slight ambiguity between the body axes and the stability 
axes, caused by the assumptions employed in thin-airfoil theory, should 
not be allowed to cause any doubt about the direction of the normal 
force. This force acts perpendicular to the plate and in a direction 
parallel to the axis, not the t axis. 

The order of presentation of the plan— form sectors in Appendix B 
of this report is a duplicate of the arrangements of the Appendixes B 
of references 1, 2, and 3. A sketch of the load distribution is pre- 
sented for each sector in Appendix B in order to provide a convenient 
visual interpretation of the load distribution. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif., Oct. 10, 19^9. 
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APPENDIX A 

SUMMARY OF EXPRESSIONS FOR LOAD DISTRIBUTION ON WING 
ELEMENTS IN TERMS OF WIND^AXES NOTATION 

Expressions Apply to Crosshatched Plan— Form Areas 

1. Supersonic leading edge in conjunction with, a subsonic leading 

edge: 
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2. Supersonic leading edge in conjunction with a subsonic trailing edge: 
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where 

E is the complete elliptic integral of the second kind with modulus 
J 1-G 2 

, = l-^ 2 nmo- J ( l-^ 2 m g ) ( l-B 2 n Q 2 ) 

B(ni— no) 

Subsonic leading edge in conjunction with a subsonic trailing edge: 



P = 


- / - 2 -- 

it v 1— Bno 


when n Q > 0; 


ni 



P = 


Bni — Gi 


Gi + k* 


(Gi-Bn 0 )(l-Gi 2 ) 1 G ^T L 


E F(q>,k) - K E (<P,k) 
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where 


Gi 


1+B 2 nin 0 — ,/( l-^ 2 ni 2 ) ( l-B 2 ^ 2 ) 
B(ni+n 0 ) 


jj» _ Gi~Bno K = F( 
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k =J 1-k' 2 E = E( 
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sin -1 


B (^° + x) 


APPENDIX B 

SUMMARY OF EXPRESSIONS FOR LOAD DISTRIBUTION IN 
TERMS OF BODY-AXES NOTATION 

Expressions apply to crosshatched and heavily 
shaded plan— form areas 


1. Triangular Wings: 
A. 



where 

E is the complete elliptic integral of the second kind with modulus 


J 1-G 2 
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q _ ( l-m 2 tan 2 3 ) + B 2 (m 2 — tan 2 3) _ 

2Bm( l+tan 2 3 ) 

«/[(l-m tan P) 2 -jB 2 (m+ tan P) g ][(l4-m tan p) 2 -B g (m— tan 3) 2 ] 

2Bm(l+ tan 2 P) 


B. 




AP\ 

q K 


when tan P = m 


tox tap 131 

P a/ ( id — “j-) ( 1+m tan 3 ) 


Pi = r 


GiBm 


EV (l-m 2 ) 


when tan 3 = 


1-Bm 

B+m 


Pi = 




l+m tan 3 


« v (l+m tan 3)+B(ni— tan 3) 
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when m <tan 3 < 


1-Bm 

B+m 


Pi = 


' [B(m+ tan 3)-Gi(l-m tan 3)](l+m tan 3) 

[B(m- tan 3)+Gi(l+m tan 3)](l-ni tan 3)(l-Gi 2 ) 


Gi+k* 


— + E + Ll V 1-0 ! 2 _ p (<p,,k)-K E (cp,k)] 

Gl *J Gx2-k»2 


where 


Gi = 


(l-ffl 2 tan 2 3)-£ 2 (m 2 — tan 2 3) 
2B( 1+m 2 ) tan 3 


J\ ( 1+m tan 3) 2 -B 2 (m-tan 3) g ][(l-m tan 3) 2 -B 2 (m+tan 3) 2 ] 

2B(l+m 2 ) tan 3 


= a/ 1— k*2 


_ Gx(l+m tan 3)+B(m— tan 3 ) 
(l+m tan 3)+GiB(m-tan 3) 


<P 


sin - 1 


Gi 2 -4c f2 

~Gik 


K = F(-,k); E = E(-,k) 




14 


NACA TN 2007 


AP 


l<x(m+tan P) 


‘l 'Ci •J B 2 (m+tan P) 2 -(l-m tan P) 2 
/ AP \ 8a(m+tan P ) 

^ ^ 'C 2 n J B 2 (m+tan P) 2 -(l-m tan P) 2 


[B(m+tan p)— (l-m tan P) ](m-tan P) 
[B(m-tan P)+(l+m tan p)](m+tan P) 


2 m[(l-B tan P)- (B+tan P) ] 

—————— ^ y 

[m(B+tan P)-<l-B tan P)](m+ ^- 7 ) 


tan - 1 
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[m(B+tan (3)— ( 1— B "bail p)](m+ ^,) 
2m[(l-B tan p)- ^ (B+tan p)] 


D. 


/ 
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/ 
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AP 


4x(m+ tan P) 


q 'l>i J B 2 (m+tan P) 2 -(l-m tan P ) 2 
8 a(m+tan p) 


(ft- 


tan' 


^ [ m(B+tan p)— (l-B tan P)](m+ -|y) 


^ 2 1 ij B 2 (m+tan p) 2 -(l-m tan P ) 2 


2 m [(l-B tan P)— (B+ tan P)] 
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E. tan 0 

B+m 



AP 


4 a(m— tan 0 ) 


^ Ei v B (m-tan 0 ) — ( l+ni t an 0 ) 2 


(fi- 

E 2 Ti 


*kx(m+ tan p) 


\/;B 2 (m+tan P) 2 -(l~ci tan ( 3) 2 


it 
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sin' 


— i 
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E*. tan 3 = — — (See footnote a.) 
B+m 
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4a(m 2 +2mB— l) 


AP\ = 

^ ip* J**/ B 2 (m 2 +2mB— 1 ) 2 — ( B +2m-m 2 B ) 


— + 


sin - 


tat 


-i 


B 2 (m 2 -t- PmB — ' 1 ) [(mB— 1 )+ (B+m) ]— (B+Pm-fl^B) [ (B+m)— (Bm— l) ] 


B(m- Ip) [(B+m) 2 +(mB-l) 2 ] 
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' ( m+2B-rB 2 m) — Ip (B 2 +2mB-l) 
(B 2 +l)(m+ pr) 


ft). - 


ha.( m 2 +PmB— l) 


E 3 * J B 2 ( m 2 +2xdB— 1 ) — ( B+2m-m 2 B ) 


a Left leading edge hits Mach cone from apex. 
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2. Swept-back wing components: 
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J* . tan (3 


Bin — 1 
B+m 


( See footnote b . ) 




^Left leading edge hits Mach cone from apex. 
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-p_ 1 

K*. tan 0 = (See footnote b.) 

B+m 



^Left leading edge hits Mach cone from apex 
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where rj* in should be replaced by tj ' + — mc r ) 

V in J be replaced by rj' - - mc r ) 

^ \ 2 
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should he replaced by i) r + (jr — ^r) 
should be replaced by t)' — (^- — mc r ) 



S. 





30 


NACA TN 2007 



b 
2 

S**. 


in 


<1 Jr 


should, be replaced by t\' 




(f )*- -(f). 

(?)„_■(?} <?}<?} 

where T| ' in f J should be replaced by tj r + — 
V q / N 2 

>■ - (fl should he replaced by r\ * — 


NACA TN 2007 


31 


T. 



U. 






4a(m— tan (B ) 

J B 2 (m— tan 3) 2 — (l+m tan 3) 2 



32 


NACA TN 2007 


V. 



AP \ _ t<x( m+tan |3) 

^ 'V J B 2 (m+taja p) 2 -(l-in tan (3) 2 
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Figure / -The triangular i trapezoidal } and rectangular pian- 
form types investigated . 
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Figure 2- Swept-back plan forms and Mach cone 
configurations investigated. 
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x > y, z rectangular coordinates of wind axes 

rectangular coordinates of stability axes 
£>y>£ r ectanguia r coordinates of body axes 



Figure 3. -Coordinate axes systems used in analysis 
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